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We present a detailed study in the paramagnetic regime of the frustrated s — 1/2 spin-compound 
linarite, PbCuS04(OH)2, with competing ferromagnetic nearest-neighbor and antiferromagnetic 
next-nearest-neighbor exchange interactions. Our data reveal highly anisotropic values for the sat- 
uration field along the crystallographic main directions, with ~ 7.6, ~ 10.5 and ~ 8.5 T for the 
a, b, and c axes, respectively. In the paramagnetic regime, this behavior is explained mainly by 
the anisotropy of the g-factor but leaving room for an easy-axis exchange anisotropy. Within the 
isotropic J1-J2 spin model our experimental data are described by various theoretical approaches 
yielding values for the exchange interactions Ji ~ -100 K and J2 ~ 36 K. These main intrachain 
exchange integrals are significantly larger as compared to the values derived in two previous studies 
in the literature and shift the frustration ratio a = J2/I Ji| f» 0.36 of linarite closer to the ID critical 
point at 0.25. Electron spin resonance (ESR) and nuclear magnetic resonance (NMR) measurements 
further prove that the static susceptibility is dominated by the intrinsic spin susceptibility. The 
Knight shift as well as the broadening of the linewidth in ESR and NMR at elevated temperatures 
indicate a highly frustrated system with the onset of magnetic correlations far above the magnetic 
ordering temperature Tn = 2.75(5) K, in agreement with the calculated exchange constants. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.30.Et, 75.30.Gw, 75.30.Kz, 75.40.Cx, 75.50.Ec, 76.30.Fc, 76.60.-k 



I. INTRODUCTION 



One-dimensional (ID) quantum magnetism involves 
the study of materials with magnetic ions of low-spin 
state, s = 1/2 or 1, which are coupled magnetically along 
one crystallographic direction only. Due to the rich na- 
ture of the low-temperature and field-induced phases in 
these materials, quasi-lD quantum magnets (Q1DQM) 
and their effective ID models have attracted much inter- 
est in the last decades from experimentalists and theo- 
rists alike. 1-4 Various Q1DQM exhibit highly fascinating 
field-induced phenomena among which the Bose-Einstein 
condensation (BEC) of one-magnons (triplons) 5-8 or two- 
magnon bound states near the saturation field 9,10 are 
particularly relevant examples. Recent attention has 
been focussed on the latter phenomenon which is well 
established for ID models with competing interactions, 
where via the inclusion of antiferromagnetic ( AFM) next- 
nearest-neighbor (NNN) interactions magnetic frustra- 
tion plays a decisive role. A prototype of these systems is 
the spin s — 1/2 chain with ferromagnetic (FM) nearest- 
neighbor (NN) and AFM-NNN interactions, 



H = Ji ^S, • S l+1 + J 2 S ; • S i+2 - h^Sf, (1) 
l 1 1 

with the NN exchange J\ < 0, the NNN interaction J2 
> 0, and h as the external magnetic field along the 
z direction. For these materials, it has been shown 
that the ground state of the ID Hamiltonian (1) has an 
instability towards field-induced multipolar Tomonaga- 
Luttingcr-liquid phases, 11-18 and which is interpreted 
as a hard-core Bose gas of multimagnon bound states 
undergoing Bose-Einstein condensation at fields slightly 
lower than the saturation field, H s . The ground state of 
the ID Hamiltonian (1) is fcrromagnetically ordered for 
a = |J 2 M| < 0.25. At \J 2 /Ji\ = 0.25, the FM state 
is degenerated with a singlet state, while for IJ2/J1I > 
0.25 the ground state is an incommensurate singlet. For 
a corresponding quasi- ID system with small but finite 
interchain couplings, the latter situation might possibly 
result in a helical magnetic structure with an acute pitch 
(screw) angle. However, depending on the strength of the 
exchange anisotropy and/or the spin-phonon coupling, 
other types of helicoidal magnetically ordered states can 
be realized, such as helical structures with obtuse pitch 



angles or cycloidal helices, ordinary Neel states, a spin- 
Peierls phase, or various other massive phases. 19 Despite 
its initially simplistic appearance, we would like to stress 
that the J1-J2 model under consideration has not yet 
been fully investigated theoretically, especially with re- 
spect to measurable physical quantities. (To illustrate 
this point see e.g. Rcfs. [20,21] and references therein.) 

The list of materials, which have been thought 
to be described well by Eq. (1) and which could 
serve as a testing ground, includes several quasi-one- 
dimensional (Q1D) edge-sharing chain cupratcs, such 
as Rb 2 Cu2Mo 3 Oi2, LiCuV0 4 , NaCu 2 2 , LiCu 2 2 , and 
PbCuS0 4 (OH) 2 . 10 ' 22 " 33 For most of these systems the 
magnetic couplings - and thus H s — are rather large pro- 
vided the system is not in the vicinity of a ferromagnetic 
critical point such as Li 2 ZrCuC>4. In the latter case, sin- 
gle crystals are still not available making it generally dif- 
ficult to study effects close to H s experimentally using 
stationary fields. 34-37 

Possibly, only one of these materials, the natural min- 
eral linarite, PbCuS04(OH) 2 , reflects an optimum com- 
promise of a comparatively low a-value not far from the 
critical point, viz., a low saturation field, and the avail- 
ability of single crystals. Despite this, linarite has only 
been studied superficially to date. According to these ini- 
tial studies, it has been proposed to represent a frustrated 
quasi-lD J\ — J 2 model magnet described by Eq. (1). 
Linarite crystallizes in a monoclinic lattice (space group 
P2 1 /m; a = 9.682 A, b = 5.646 A, c = 4.683 A, /3 = 
102. 66°), 38 in which Cu0 2 units are aligned in chain- 
like structures along the b direction (Fig. 1). Each of the 
Cu 2+ atoms is fourfold coordinated, with the surrounding 
oxygens forming a flat tetragon close to a square, the cele- 
brated CuC>4 6 plaquette-" brick" of all undoped cuprates. 
This is supplemented by four hydrogen atoms as ligands 
to these planar oxygens. The Cu coordination is com- 
pleted by two further oxygen atoms from the SO4 groups, 
yielding a distorted octahedron and a slightly non-planar 
"waving" (buckled) Cu0 2 -chain structure. This is in con- 
trast to the much better studied planar counterpart of the 
edge-shared chain cuprates LiCu 2 2 and LiVCuCv 

Initial susceptibility and zero-field specific heat data 
have been interpreted in terms of a dominant magnetic 
coupling along the chain, with a predominant FM-NN 
interaction J x = -30 K and a weaker AFM-NNN cou- 
pling J 2 = 15 K (a = 0.5). 33 This way, a strong coupling 
scenario in terms of interacting and interpenetrating sim- 
ple AFM Hciscnbcrg chains would be envisaged for the 
J1-J2 model under consideration. However, more recent 
studies by Yasui et al. 22 indicate a rather different weak- 
coupling regime, namely: J\ = (-13 ± 3 ) K and J 2 = (21 
± 5)K (a = 1.6), which have been obtained from a fit 
to the susceptibility data using a high-temperature ex- 
pansion up to the fourth-order in the temperature range 
50 < T < 350 K. Such a weak-coupling regime should 
give rise to more pronounced frustration and fluctuation 
effects. 39 The observation of a magnetically ordered state 
below Tn ~ 2.8 K has been discussed in terms of a possi- 




FIG. 1: (Color online) The crystallographic structure and the 
main exchange paths in PbCuS04(OH)2- The Cu02 units are 
aligned in the be plane, forming edge-sharing CuC>2 chains 
along the b direction. In order to illustrate the coordination of 
the S atoms, oxygen tetrahedra are highlighted in the sketch. 
The two inequivalent alternating "left" and "right" proton 
positions distinguished by bond lengths and bond angles are 
according to Ref. 38. 



ble helical ground state with an acute pitch angle. 33 The 
helical nature of the ground state is supported by a re- 
cent study of the dielectric constant in this material. 22 
However, microscopic studies such as NMR or neutron 
scattering measurements to prove these predictions are 
still lacking. Although in both references a similar con- 
clusion has been obtained regarding the basic nature of 
the magnetic ground state of linarite, 22,33 the physical 
characteristics of the spiral and thus the magnetic prop- 
erties would differ enormously due to e.g. the relevance of 
quantum effects on the pitch angle and on the magnetic 
moment in the ordered state. 

In view of these conflicting results, we started a de- 
tailed study on this material combining macroscopic and 
microscopic experimental techniques with different theo- 
retical methods to resolve the magnetic exchange param- 
eters of PbCuS04(OH) 2 and, thus, its underlying ground 
state. Measurements of the static susceptibility and the 
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saturation magnetization have been performed. The in- 
trinsic spin susceptibility is investigated utilizing ESR 
and NMR. Measurements of the g-factor as well as of 
the ESR and NMR linewidths indicate an appreciable 
magnetic anisotropy in our system and further empha- 
size the highly frustrated character of linarite with an 
onset of magnetic correlations far above the magnetic 
ordering temperature. An analysis of our data within 
advanced theoretical methods, such as density matrix 
renormalization group (DMRG), hard-core boson tech- 
nique and local (spin-polarized) density approximation 
(L(S)DA+J7) calculations, which all take into account 
the nontrivial interplay of quantum effects and frustra- 
tion beyond linear spin- wave theory, yield new values for 
the exchange couplings along the chain directions of J\ 
~ -100 K and J2 ~ 36 K. These values are substantially 
larger than those determined previously. 22,33 Within our 
extensive analysis, we succeeded in estimating the order 
of magnitude of the interchain couplings. The theoreti- 
cal part of the present paper should be understood also 
as a step in the direction to fill the existing gap of the 
not fully explored J\ — J2 model, beyond providing an 
assignment of the main exchange parameters for a spe- 
cific material only. In particular, a precise knowledge of 
the main exchange parameters is one of the prerequisites 
to attack in a realistic manner rather complex and not 
yet fully understood phenomena such as multiferroicity 22 
and multipolar phases/spin nematics 12 reported or sug- 
gested for the title compound, too. 



II. METHODS 
A. Experimental 

1. Samples and diffractions 

The single crystals of PbCuS04(OH)2 used in this 
study for the magnetization, NMR, and X-band ESR 
measurements are natural minerals with their origin in 
California, USA (origin 1: Blue Bell Mine, Baker, San 
Bernadino). A second set of naturally grown single crys- 
tals of smaller size (origin 2: Siegerland, Germany) has 
been used for the ESR measurements in a resonant cav- 
ity at a frequency of about 93 GHz. All crystals show 
well-defined facets and the principal axes b and c can be 
identified easily. Single crystallinity of our samples has 
been checked by x-ray diffraction. For both sets of sin- 
gle crystals no magnetic impurity phases were observed 
within experimental resolution, as evidenced by the ab- 
sence of a low-temperature Curie tail in the magnetic 
susceptibility. For all measurements the samples were 
oriented along the principle crystallographic directions 
with a possible misalignment of less than 5°. 



2. Magnetization 

Temperature-dependent magnetization studies of 
linarite were performed using a commercial SQUID mag- 
netometer in the temperature range 1.8 - 400 K and in 
an external magnetic field of 0.4 T. Magnetization curves, 
M(H), were measured at T = 1.8 and 2.8 K in a com- 
mercial Physical Properties Measurement System using a 
Vibrating Sample Magnetometer (VSM). Magnetization 
data were collected while sweeping the magnetic field 
using sweep rates of about 50Oe/s for both increasing 
and decreasing field regimes. Note that due to hystere- 
sis around the phase transitions observed for the 1.8 K 
data, the sweep rate was significantly varied in these re- 
gions in order to check for sweep-rate dependent effects 
on the phase transitions. Using quasi-static conditions, 
the observed small hysteresis in the M(H) curves became 
negligible, as shown below. 



3. Electron Spin Resonance (ESR) 

For the ESR experiments two different setups were 
used. The commercial X-band spectrometer operates 
at a frequency of about 9.6 GHz. This allows sweeps of 
the magnetic field up to 1 T. The setup is equipped with 
a continuous- flow liquid- helium cryostat, enabling mea- 
surements from room temperature down to about 4K. 
The cryostat is inserted in a rectangular microwave res- 
onator in the TE102 mode configuration. Samples are 
mounted on a quartz sample holder, which is centered 
in the resonator at the maximum of the microwave mag- 
netic field. The sample holder can be rotated with re- 
spect to the external magnetic field. By using an ad- 
ditional external magnetic modulation field, the lock-in 
detection technique is applied. The second setup consists 
of a microwave vector network analyzer (MVNA) and a 
15 Tesla superconducting magnet. 40 This setup allows 
phase-sensitive measurements at different frequencies in 
the range from 30-800 GHz. Most experiments were car- 
ried out using a home-built cylindrical resonator in the 
TE011 configuration with a resonance frequency of about 
93 GHz. For this, the samples were mounted on a quartz 
needle in the center of the resonator. The resonator is 
coupled to a metallic waveguide, which is placed in the 
center of the superconducting magnet. The very high 
quality factor of the order of Q ~ 10 4 - and, thus, the 
high sensitivity of this resonator setup - allowed the mea- 
surement of single crystals of linarite with dimensions of 
approximately 0.5 x 0.5 x 1 mm 3 with very high signal- 
to-noise ratio. 



4- Nuclear Magnetic Resonance (NMR) 

1 H-NMR ( 4 7 = 42.5749 MHz/T) and 207 Pb-NMR 
( 207 7 = 8.90 74 MHz/T) measurements were performed 
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using a phase-coherent Tecmag spectrometer with a He- 
flow cryostat for temperatures down to 4.2 K and mag- 
netic fields of 2 and 4 T, respectively. Temperatures be- 
low 4.2 K were achieved by pumping on the helium bath. 
The NMR spectra were determined using a 7r/2 - r - tt 
Hahn spin-echo pulse sequence. Special care was taken 
to avoid extrinsic signals from parasitic 1 H atoms around 
the sample. Due to the strong increase of the linewidth 
at low temperatures, additional field sweeps at constant 
frequency have been performed. This way, we can ensure 
that neither the selected frequency excitation spectrum 
by the pulse width (typically a 7r/2-pulse is about 3 fj,s) 
nor the quality factor of our coil do artificially narrow 
the detected lines of PbCuS04(OH)2. The spin-lattice 
relaxation time Ti has been recorded using an inversion- 
recovery pulse sequence (it - r var - %/2 - r - 71") with 
variable delay r var and a Hahn spin-echo detection at 
the end. Typical conditions of excitation were 3 /is and 
6 /is for a ir/2- and 7r-pulse, respectively. Repetition rates 
were in the range 100 - 400 ms despite short spin-lattice 
relaxation times T\ in order to avoid any local heating 
at the sample site. For H \\ b, the spin-lattice relaxation 
rate of the 1 H nucleus was determined using a saturation- 
recovery sequence with an echo subsequence at the end, 
i.e., (7r/2 - Tdei) n - T var - 7r/2 - t - 7t with the delay time 
Tdei and n as the number of repetitions of the first cycle. 
Calibration of the fields has been performed using the 
1 H- and 2 D-NMR resonance frequencies of hydrogenated 
and dcutcratcd water at room temperature for the 2 and 
4T experiments, respectively 

B. Theoretical methods 

1. DMRG and TMRG 

Initially, we analyzed the saturation field using well- 
known rigorous expressions valid in the so-called two- 
magnon and onc-magnon sectors depending on the 
strength of the interchain coupling. Next, we considered 
the magnetization curve at low T, say 1.8 K, for H _L 
(be) and compared our calculations with the experimental 
data shown in Fig. 3 for this direction. To calculate the 
magnetic susceptibility we employed the transfer-matrix 
renormalization group (TMRG) method. 41-43 In our cal- 
culations, 80-160 states were retained in the renormal- 
ization procedure and the truncation error was less than 
10~ 4 down to T = 0.003 1 Ji | . 

To calculate the static spin-structure factor S(q) = 

(1/i 2 ) Ey=i [<5f S/) - (Sf){S*)] , where Sf denotes the 
z-component of the spin at site i, we used the density- 
matrix renormalization group (DMRG) method. 44 We 
studied single chains (two coupled chains) with length 
up to L = 512 (L = 96) and keeping up to m = 4000 
(m = 2000) density-matrix eigenstates in the renormal- 
ization procedure such that the truncation error was less 
than 10 -9 (10 -6 ). In the single-chain case, the calcu- 
lated values were extrapolated to the thermodynamic 



limit L — > oo. We calculated the magnetization curve 
at very low temperature (T/J <C 1) using the DMRG 
technique. The magnetization for a given magnetic field 
was obtained as 

where E n and ip n are the n-th eigenenergy and the eigen- 
state (n = denotes the ground state), respectively, kg 
is the Boltzmann constant, and m c is a cutoff number 
of the excited-state energies. The cutoff number m c de- 
sirably should be sufficiently smaller than the number 
of states m kept in the density-matrix renormalization 
step. In this paper, for a fixed system length L = 64 and 
temperature T = 1.8 K, m c varied from 40 to 100 while 
keeping m = 1000 and M was extrapolated to the large 
m c limit. 

Applying the TMRG technique we analyzed the mag- 
netic spin susceptibility, x(T), to obtain its maximum 
position, T™ ax , as a function of the frustration parame- 
ter a. Notice that in the adopted random phase approx- 
imation (RPA) for the interchain coupling (IC) T™ ax of 
the 3D susceptibility, Xsd(T), is independent of the IC, 
the (/-factor, and the background susceptibility %o- 

2. Local spin density L(S)DA+U calculations 

In the second part, we calculate the total energies for 
various prepared magnetic states, i.e., a ferromagnetic, 
and various antiferromagnetic states, whose total energy 
differences were mapped onto those of corresponding spin 
states of a generalized J\-Ji model with supplemented in- 
terchain interactions. This way, the main exchange inte- 
grals could be extracted. 36,45 The density functional the- 
ory (DFT) based electronic-structure calculations were 
performed using the full potential local orbital scheme 
FPLO9.00-33. 46 ' 47 Within the scalar relativistic calcula- 
tion the exchange and correlation potential of Perdew 
and Wang was applied. 48 For the LSDA+C/ calculations 
we varied U in the physical relevant range from 5 to 8 
cV using the mean-field approximation of the double- 
counting correction. To ensure convergency we consid- 
ered 518 k points within the irreducible part of the Bril- 
louin zone. 



III. MAGNETIC SUSCEPTIBILITY 

Figure 2 displays the temperature dependence of the 
macroscopic susceptibility x(^) °f linaritc with the mag- 
netic field applied along the three principle crystallo- 
graphic directions. The susceptibility exhibits two char- 
acteristic features in the low-temperature region, i.e., a 
maximum around T™ ax = 4.9 ± 0.3 K (averaged over all 
three crystallographic directions, see below), just above 



the Neel temperature Tn visible as a pronounced kink 
around 2.8 K. While the maximum around 5K is very 
common in quasi-lD frustrated "ferromagncts" in the 
vicinity of the critical ferromagnetic-helical point 3 (see 
also Fig. 13 in Sect. VI) and is associated to low-lying fer- 
romagnetic excitations, the kink at lower temperatures is 
related to the transition into a long-range magnetically 
ordered state. From the derivative d(xT)/dT of the sus- 
ceptibility data (insets in Fig. 2) the transition temper- 
ature T N = 2.75(5) K is determined. 

While the qualitative behavior of x(T) with respect to 
the three crystallographic directions is the same, the ab- 
solute values oixmax at T™ ax differ slightly with T™ ax = 
5.0 ± 0.2 K for H || a, T™ ax = 4.6 ± 0.3 K for H \\ b, and 
r™ = 4.9 ± 0.2 K for H \ \ c. From this observation, 
one already might speculate that the 6-axis should be the 
easy axis, despite of the larger error bars of T™ ax . We 
would like to stress, that our T™ ax values are in agree- 
ment with those reported in Ref. 33, while they differ 
from those of Yasui et al., 22 where T™ ax lies between 3.8 
- 5 K and becomes maximal for if | a' . The origin of this 
discrepancy is not clear up to now, although one possible 
scenario could be related to small impurity contributions 
yielding an additional Curie-like contribution, this way 
shifting T™ ax to lower temperatures. 

In Fig. 3, we present the magnetization data M({j,qH) 
and the derivatives dM /d(iioH) of PbCuS04(OH)2 as a 
function of field at T = 1.8 K < T N and T = 2.8K > 
Tn with H _L (be) (hereafter named a±), H \\ b, and H 
|| c. 49 Although here we will focus on the properties of 
linaritc in the paramagnetic regime, we added the low- 
temperature data in order to extract the saturation field, 
H s , in the case of reduced thermal fluctuations. Mea- 
surements were made for increasing and decreasing field; 
no hysteresis has been observed within the experimen- 
tal resolution as result of a low sweep rate. Compar- 
ing the magnetization versus field (and its derivative), a 
large anisotropic response is observed with M s ^ a± w 1.16 
/l*b/Cu atom, M s ^ sa 1.05 //b/Cu atom, and M S:C rj 1.15 
/ie/Cu atom for the a±, b, and c direction, respectively. 
Together with the anisotropic values of the saturation 
field fJ, H s , i.e., [i H s , a± ~ 7.6 T, ^ #s,6 ~ 10.5 T, and 
fioH StC ps 8.5 T, this implies an appreciable anisotropy of 
the magnetic exchange in our compound. 

The anisotropic magnetic behavior of PbCuS04(OH)2 
is also evidenced in the number of transitions observed 
in the dM/d(/j,oH) curve at 1.8 K. While for the a± and 
c directions only one sharp peak indicative of a phase 
transition is observed (around fJ-oH C 3 ~ 5.2-5.5 T), for 
the b direction, i.e., for the chain direction, three clear 
transitions could be resolved at /-LqH^ ps 2.7 T, fJ-oH^ 2 w 
3.3 T, and HoH^ 3 « 5.8 T. These transitions at low fields 
possibly assign rotation and spin-flop reorientations out 
of a helical ground state for the field aligned along the 
easy magnetic axis. We will discuss this behavior in more 
detail in a forthcoming paper. 50 

Both the anisotropy and the overall small saturation 
field of linarite are very important and unique features 
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FIG. 2: (Color online) Temperature dependence of the mag- 
netic susceptibility of PbCuS04(OH)2 for an external mag- 
netic field of 0.4 T applied parallel to the a, b, and c axes. In 
the insets, the derivatives d(xr)/dT as function of tempera- 
ture are shown, from which the magnetic ordering tempera- 
ture Tn ~ 2.75(5) K can be evaluated. 



with respect to the predicted exotic high-field phases. 
This makes PbCuS04(OH)2 suitable for the investiga- 
tion of such phenomena in easily available static mag- 
netic fields. In order to resolve the origin of the magnetic 
anisotropy in linaritc, which could either stem from g- 
factor anisotropy, anisotropic exchange, or antisymmet- 
ric Dzyaloshinskii-Moriya interactions to name a few, we 
performed ESR investigations. 



IV. ELECTRON SPIN RESONANCE (ESR) 

ESR measurements were performed for different single 
crystals of linarite. A small single crystal from origin 2 
was used for the measurements in the resonant cavity at 
a frequency of about 93 GHz. At this frequency, the res- 
onance field is found around 3 T. The spectrum consists 
of a single line of Lorentzian shape (Fig. 4). From a fit to 
the ESR lines, the intensity, the resonance field, and the 
linewidth are extracted. From those parameters the inte- 
grated ESR intensity, which is determined by the intrinsic 
spin susceptibility, was calculated (sec Fig. 9 and discus- 
sion in the following NMR section). Since the resonator 
was inserted into the cryostat its quality factor Q and 
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its resonance frequency v TCS change with temperature. Q 
and u res are determined with frequency sweeps around 
the resonance frequency at every temperature and the 
spin susceptibility Xs and resonance field H res can be 
corrected accordingly. 

The temperature dependencies of the ESR resonance 
field and ESR linewidth are plotted in Fig. 5. 51 In the 
high-temperature regime (roughly above 50 K) the res- 
onance fields are temperature independent. A shift of 
the resonance position of an ESR signal as a function of 
temperature is associated with the development of inter- 
nal magnetic fields in the system. Along the b direction a 
shift is observed only at low temperatures close to the or- 
dering temperature Tn- This shows that an internal field 
develops along the chain direction only when the actual 
3D ordering occurs. However, along the a and c direc- 
tions a shift in the resonance field is observed already at 
much higher temperatures starting at around 50 K and 
developing smoothly with decreasing temperature. This 
indicates that the Cu spins arc primarily aligned in the 
ac plane, but do not point along the chain direction in 
the paramagnetic regime. The internal fields developing 
are, therefore, predominantly directed perpendicular to 
the chain direction. At first glance, the above consider- 
ations are in some conflict with the magnetization mea- 
surements indicating that the b axis might be the easy 
axis of our system. However, taking into account that 





FIG. 3: (Color online) Magnetization curve M(fioH) and its 
derivative AM /d(fj,oH) for magnetic fields applied along a±, 
b, and c of PbCuS0 4 (OH) 2 measured at T = 1.8 K < T N 
(black curves) and T — 2.8 K > T N (red curves). 
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FIG. 4: (Color online) ESR spectra for different temperatures 
measured at 93 GHz with the magnetic field H applied along 
the b and c axes. Spectra for H || a (not shown) are sim- 
ilar to spectra for H || c. Arrows on the left side indicate 
temperatures of the individual ESR spectra. 



the magnetic field used in the ESR measurements (~ 3 
Tcsla) will probably be strong enough to rotate the pre- 
spiral orientation out of the easy-plane (see Sect. Ill), 
even at high temperatures above Tn short-range ordered 
clusters could explain our observations with predominant 
internal fields perpendicular to the chain direction. 

From the resonance field, H res . the effective g-factors 
along the crystallographic directions can be determined 
as j = h ■ v res /{^H res ). For the high-temperature 
regime the effective ^-factors are found to be g a = 2.34, 
g b = 2.10, and g c = 2.28. 

Having established the (/-factors for the three princi- 
pal crystallographic directions, we can analyze the mag- 
netic anisotropy in our system. In Fig. 6, we present the 
spin expectation value < S z >= M/Ng^s (N: Avogadro 
number) of PbCuS04(OH)2 as a function of the scaled 
field g/J-oH at 2.8 and 1.8 K and for the three crystallo- 
graphic directions a±, b, and c as derived from the ex- 
perimentally determined magnetization data M(H) (Fig. 
3). The extracted spin expectation value corresponds 
to the Cu spin 1/2. In the paramagnetic regime above 
Tn, the anisotropy of both saturation magnetization and 
saturation field is explained mainly by the anisotropy 
of the <?-factor. Note that there is a difference in the 
calculated and directly measured saturation magnetiza- 
tion for H || a± due to the g- factor. For the ESR ex- 
periment the magnetic field was aligned along a, while 
for the magnetization measurement the field was aligned 
perpendicular to the be plane. For an alignment along 
a the saturation magnetization would be slightly larger 
and would match the calculated value. For temperatures 
smaller than Tn, however, the anisotropy cannot be de- 
scribed by the g-factor anisotropy. Additional contribu- 
tions from symmetric exchange anisotropy and/or possi- 
bly Dzyaloshinskii-Moriya interactions need to be taken 
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FIG. 5: (Color online) (a) Resonance fields at about 93 GHz 
and (b) linewidths for resonance fields of about 0.3 and 2.9 T 
of the ESR signals as a function of temperature for the exter- 
nal magnetic field applied along the crystallographic axes a, b, 
and c. The resonance field is temperature independent above 
~ 50 K and decreases with decreasing temperature. For the 
high-temperature resonant fields the corresponding g-factors 
are listed. The linewidth is significantly smaller for H || b. 



into account. 

As shown in Fig. 5, the ESR linewidth is strongly 
anisotropic. At a frequency of 93 GHz corresponding to 
a resonance field of about H res w 3T, the linewidth 
for H || b is with AH « 0.13 T much smaller than 
AH w 0.7 T for H \\ a and c. For all directions the 
linewidth is almost constant or only weakly dependent 
on temperature above 50 K. At the X-band frequency of 
about 9.6 GHz - corresponding to H res « 0.3 T - the tem- 
perature dependencies of the ESR signals show a similar 
behavior (squares in Fig. 5(b)). The linewidths along a 
and c are - within the uncertainty of the measurement - 
identical to the linewidths at higher fields. The decreas- 
ing intensity of the ESR signal with increasing tempera- 
ture makes it difficult to analyze the lines perfectly along 
those orientations up to room temperature, however, the 
linewidth appears to stay constant. For the b direction 
at high temperatures, the linewidth is constant and fairly 
narrow with AH w 0.25 T, which is unexpected since it is 
broader than for the larger field of about 3 T. This effect 
can be explained by the strongly anisotropic linewidth to- 
gether with a slight misalignment of the sample of about 
5° for this particular measurement. 

On approaching temperatures below 50 K a broaden- 
ing of the lines is observed for both fields. The small 
resonance field of about 0.3 T limits the reliability of the 
measurements for fields aligned along a and c, since the 
linewidth exceeds the resonance field and a fit to the data 



cannot be done accurately anymore. However, the values 
are close to the linewidths for the ten times larger field 
of about 3 T. 

The ESR linewidth depends on (dipolar) spin-spin or 
anisotropic exchange interactions as well as on the de- 
velopment of internal magnetic fields. The fact that the 
linewidth is the same for different fields indicates that 
inhomogeneous broadening effects arc rather small. As 
we approach lower temperatures the spin-spin correla- 
tion length increases and short-range magnetic correla- 
tions develop. The change of linewidth as a function of 
temperature yields information about the dimensionality 
and type of interactions in the system. The broadening 
can be analyzed in terms of 



AH(T) = AH + C[(T - T N )/T N ]- p 



(3) 



as the ordering temperature is approached. The 
linewidth, AH(T), is divided into a non-critical constant 
part, AHq, and a temperature-dependent critical part, 
AH cr i t (T). 52 The exponent p yields information about 
the effective dimensionality of the correlated spin sys- 
tem and its change by approaching a long-range-ordered 
ground state. Fits to our data over the low-temperature 
range with a fixed Tn = 2.0 K at 3T result in critical 
exponents p = 0.5-0.8 (Fig. 7). Note, that the magnetic 
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FIG. 6: (Color online) The experimentally determined spin 
expectation value < S z > = M/Ngfi B of PbCuS0 4 (OH)2 as 
a function of the scaled magnetic field g^oH along different 
crystallographic directions and temperatures (a) T = 2.8 K 
> T N and (b) T = 1.8 K < T N . 
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ordering temperature Tn = 2.0 K was chosen according 
to the field dependence of T N as determined by ther- 
modynamic bulk measurements. 50 For a ID Hciscnberg 
magnet the critical exponent was found to be about p 
= 2. 5. 53 However, on approaching the ordering temper- 
ature this value can change significantly and values of p 
ss 0.6 have been reported. 52 Such a critical exponent is 
interpreted as signalling the appearance of 3D antiferro- 
magnetic fluctuations. The fact that the broadening in 
our system occurs already at temperatures around 50 K, 
i.e., 15 times higher than the actual ordering tempera- 
ture, points to appreciable magnetic fluctuations at el- 
evated temperatures indicative of substantial interlayer 
correlations well above Tn- The actual 3D ordering at 
much lower temperatures then indicates a strongly frus- 
trated system with competing interactions on the energy 
scale of ~ 50 K. 



V. NUCLEAR MAGNETIC RESONANCE 
(NMR) 

Our NMR experiments were conducted using 1 H- and 
207 Pb as probing nuclei. Both the hydrogen and the lead 
ions occupy low-symmetric crystallographic sites with re- 
spect to the magnetic Cu sites, i.e., none of the ions 
are located exactly between two Cu ions neither along 
the chain nor between two neighboring chain structures. 
While there is only one crystallographic site for Pb, two 
incquivalent H sites can be expected according to recent 
structural investigations via neutron diffraction empha- 
sizing two different kinds of hydrogen bondings in linar- 
ite. 38 Henceforth, a single 207 Pb and two 1 H-NMR lines 
can generally be expected in our experiment. We name 
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FIG. 7: (Color online) Double logarithmic plot of the critical 
linewidth part AH cr u (from ESR) as a function of tempera- 
ture for H parallel to the three crystallographic directions a, 
b, and c. Fits to the data for temperatures up to 25 K are 
shown. 



the hydrogen site with the stronger bonding Hi, i.e., 
probably the one with the larger hyperfine coupling, and 
the one with the weaker bonding H2. For magnetic fields 
H || b, however, two nearly overlapping *H resonance 
lines imply that both hyperfine couplings are nearly iden- 
tical for this field direction. We have also tried to detect 
the 63,65 Cu spin-echo signals, but did not succeed. We 
attribute the lack of copper signal to very short spin-spin 
relaxation times T2 of linarite, which are of the order of 
20 ^s at room temperature even at the 207 Pb and *H 
sites. 

The interest in probing both *H and 207 Pb lies in the 
different coupling of the two nuclei with their neigboring 
atoms and thus different distances and symmetries with 
respect to the magnetic Cu 2+ ions. From the crystal- 
lographic structure and chemical-bonding scheme it can 
be expected that due to the distance between Pb and 
neighboring magnetic Cu ions the hyperfine coupling at 
the 207 Pb site will be dominated by dipolar couplings 
between Pb nuclei and Cu spins. At the 207 Pb site, 
the dipolar fields will be predominantly given by the 
spins of the two nearest magnetic ions, i.e., two Cu spins 
along the chain (6) direction. But also couplings to next- 
nearest neighboring Cu ions along the chain as well as 
neighboring chains along c and a can be expected to re- 
sult in small additional contributions. On the other hand, 
at the 1 H site, the effective local fields at the probing 
nuclei are probably composed of the dipole fields of sur- 
rounding magnetic Cu 2+ moments and of so-called con- 
tact fields. The latter are due to the direct neighboring 
environment of hydrogen and oxygen atoms, which me- 
diate the magnetic superexchange between magnetic Cu 
ions. In this situation a small polarization of the hydro- 
gen atoms can be envisaged. Since the hydrogen atoms 
are located very close to the be plane of Cu 2+ ions, the 
hyperfine fields will be predominantly given by the spins 
of the four nearest magnetic ions of the nearest neighbor- 
ing be plane, i.e., two neighboring Cu spins from one and 
two neighboring Cu spins from a second Cu-chain shifted 
by the lattice constant c. 

The spin-echo signal of the 1 H and 207 Pb lines was 
observed in the temperature range 5 - 400 K. We first 
analyze the resonance shift of the : H and 207 Pb-NMR 
lines for an applied magnetic field of 2 and 4T, respec- 
tively (Fig. 8). The different field values have been used 
in order to obtain reasonable frequency ranges for both 
nuclei. Due to strong transverse magnetic short-range 
correlations leading to a very short spin-spin relaxation 
time T2 of less than 5 /xs at low temperatures, a wipe-out 
of the 207 Pb-NMR signal occurs below - 10 K. The NMR 
shift is defined as the normalized difference between the 
observed resonance frequency, u> res , and the calculated 
value for the bare nucleus, 



K(T) = 



Wres — J^qHq 



(4) 



7 Mo Ho 

7 being the gyromagnetic ratio of the nucleus and Hq 
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FIG. 8: (Color online) The 1 R- and 207 Pb-NMR shift of 
PbCuS04(OH)2 in the temperature range 5 - 400 K for all 
three crystallographic directions. While the 1 H data have 
been determined in an external field of 2 T, the 207 Pb-NMR 
shift has been measured in an applied field of 4T. In the 
insets, the NMR shifts are plotted as function of the macro- 
scopic bulk susceptibility. The lines represent linear fits to 
the experimental data; for details see text. 



being determined from the 1 H and 2 D-NMR resonance 
frequency of (dcutcratcd) water at room temperature. 
From Fig. 8 a strong paramagnetic increase of both the 
1 H- and 207 Pb-NMR shift is observed with decreasing 
temperature, which arises from the interactions between 
the probing nuclei and the surrounding electrons. 

Generally, the NMR shift K to t(T) can be divided into 
two contributions, 



Ktot(T) — K sp i n (T) + K 



L orbi 



(5) 



where K sp i n {T) = Ax S pin{T) arises via a hyperfine cou- 
pling to the electronic spins and K orb stems from a 
temperature-independent orbital magnetization induced 
at the nucleus site. Here, A is the hyperfine coupling con- 
stant, which can either have a positive or negative sign, 
leading to positive or negative temperature dependencies 
of the NMR shift. From this equation it is already obvi- 
ous that NMR has an advantage over bulk susceptibility 
investigations. Via NMR one accurately measures the 
intrinsic spin susceptibility, Xspin{T), without suffering 
from temperature-independent diamagnetic core or Van 
Vlcck contributions, from free spins (impurities) and ex- 



trinsic foreign phases, which limits the accuracy of bulk 
susceptibility measurements. Therefore, it is more reli- 
able to extract the magnetic parameters from the tem- 
perature dependence of the NMR shift rather than from 
bulk susceptibility. The conventional scheme of correlat- 
ing the NMR shift K to t{T) and the bulk susceptibility 
x(T) is to plot both parameters as K to t(x) with temper- 
ature being an implicit parameter. This way, the slope 
yields the hyperfine coupling constant A, while K orb re- 
sults from the intersect at \ = 0. 

The NMR shift as function of the bulk susceptibility 
of linarite is shown in the insets of Fig. 8 for both nu- 
clei and for the magnetic field applied along the three 
crystallographic directions. Clearly, both physical prop- 
erties scale with each other for all cases in the full tem- 
perature regime. A linear fit to this data yields highly 
anisotropic hyperfine coupling constants for both 1 Id- 
nuclei, i.e., A Hl _ H \\ a = -1.2 kOe/yLtB and A H2 , H \\ a = - 
2.6 kOe/^B, A HuH \\ b = -0.9 kOc/^ B , A H ^ H \\ b = -0.7 
kOe//i B , A HuH \\ c = 1.9 kOc/yUB and A H2M]]c = -0.8 
kOe//iB. The different hyperfine couplings of the two 
inequivalent H atoms strongly support the notion of dif- 
ferent hydrogen bondings to neighboring oxygen sites in 
linarite as determined previously. 38 In contrast to these 
anisotropic values for the 1 H- nuclei, the 207 Pb hyperfine 
couplings arc dominated by a large positive isotropic con- 
tribution, which is complemented by a small anisotropic 
(dipolar) component for the three different axes, yield- 
ing overall large and positive values of Ap b u\\ a = 26.5 
kOe/// B , A PbtH \\ b = 17.9 kOe/jUs, and A PbM \\ c = 21.3 
kOe//iB for fields aligned along a, b, and c, respectively. 

After having determined the hyperfine couplings for 
the two different nuclei, the intrinsic spin susceptibil- 
ity of PbCuS04(OH)2 can be evaluated via Xspin(T) = 
(K tot (T) — K orb )/A for the three crystallographic direc- 
tions. Fig. 9 depicts these physical properties plotted 
as the inverse spin susceptibility, xJjAn' as a function of 
temperature as derived from the 207 Pb-NMR data. 54 For 
a comparison we added the spin susceptibility obtained 
from our ESR investigations after normalizing these data 
to the value of the static susceptibility at 300 K. From this 
figure it is clearly visible that the intrinsic spin suscepti- 
bility from ESR and NMR scales nicely and is practically 
identical over the whole temperature range for all crys- 
tallographic directions. 

Then, from a linear fit of the inverse susceptibility to 
a Curie- Weiss law X^p in {T) oc (T - 6cw) m the T region 
250 - 400 K, the Curie- Weiss temperature, 0cw- could be 
determined for the three directions. Its absolute value is 
isotropic within the experimental error bars as expected 



for a Cu 



2 + 



1/2 system, yielding 6 C w = 27(2) K. 



The positive value of the Curie- Weiss temperature in- 
dicates the predominance of a ferromagnetic coupling. 
Comparing the Curie- Weiss temperature to the ordering 
temperature, a ratio Ocw/Tk ~ 10 is extracted for linar- 
ite. This quantity is commonly used to judge the level 
of frustration in a compound, since frustration tends to 
suppress long-range order. A ratio of 10 shows that frus- 
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tration is definitely an important issue that needs to be 
considered in this compound. 

In Figs. 10 and 11, the linewidths of the 207 Pb- and the 
1 H-NMR spectra, respectively, are shown as function of 
temperature for all three crystallographic directions a, b, 
and c. In the insets, some 207 Pb and 1 H-NMR spectra are 
depicted for applied magnetic fields of 4 and 2 T, respec- 
tively, and for 3 different temperatures between 10 and 
150 K. In the paramagnetic phase, the spin-echo NMR 
signal has a rather isotropic Lorentzian lineshape. Since 
the linewidth at high temperatures is very small, i.e., 
about 10-25 kHz for all spectra and for all three crystal- 
lographic directions, we suggest that our single crystal is 
of rather high-quality. Note that also for the 1 H-NMR 
spectra for H \ \ b two Lorentzian lines have been used to 
fit the data since the two NMR lines do not perfectly over- 
lap at low temperatures. In Fig. 11, however, the average 
of the linewidth of both 1 H-spectra has been plotted for 
H\\b. 

Similar to the ESR linewidth, the temperature depen- 
dence of the full width at half maximum (FWHM) of the 
NMR spectra is expected to give access to the dynam- 
ics of the magnetic correlations and thus to the dynam- 
ical critical properties in the paramagnetic regime upon 
approaching Tn. The NMR linewidth is related to the 
spin-spin relaxation time, T2, and thus probes the trans- 
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FIG. 9: (Color online) Inverse spin susceptibility of linar- 
ite, X^inCT), for the external magnetic field applied along a, 
b, and c as determined via 207 Pb-NMR and ESR. The spin 
susceptibility as determined via ESR at about 3T was nor- 
malized to the high-temperature (300 K) value of the static 
susceptibility. The lines represent linear fits to a Curie- Weiss 
law in the high-temperature range from 250 to 400 K. 



FIG. 10: (Color online) The 207 Pb-NMR linewidth of 
PbCuS04(OH)2 as a function of temperature in an external 
magnetic field of 4 T applied along the three crystallographic 
axes o, b, and c. The lines are guides to the eyes. The insets 
show representative 207 Pb-NMR spectra for different temper- 
atures. The different intensities at these temperatures are not 
to scale; for details see text. 



verse component of the two-spin correlation function and 
the temporal spin fluctuations of the magnetic system 
near the critical temperature. Considering a compound 
with anisotropic magnetic exchange interactions such as 
linarite, it can be expected that the NMR linewidth is 
dominated by spin fluctuations along the magnetic easy 
axis, with spin fluctuations perpendicular to the easy axis 
only contributing to the non-critical broadening. Hence, 
taking into account that the linewidth probes transverse 
spin fluctuations the broadening of the NMR line should 
be most prominent for magnetic fields perpendicular to 
this (easy) axis. 

For all NMR spectra a pronounced broadening of the 
line has been observed below ~ 75 K for the 207 Pb signal 
and below ~ 50 K for the 1 H spectra. This broadening 
points to short-range correlations developing already at 
temperatures T 3> Tn- Comparing both the response at 
the two different nuclei and for the three different crystal- 
lographic directions, one can easily see that (i) the broad- 
ening of the NMR line is shifted to lower temperatures 
but slightly enhanced for the X H spectra and (ii) that 
particularly for the 2T X H-NMR data the broadening is 
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FIG. 11: (Color online) The 1 H-NMR linewidth of 
PbCuS04(OH)2 as function of temperature in an external 
magnetic field of 2 T parallel to the three crystallographic axes 
a, b, and c. The lines represent guides to the eyes. The insets 
show representative 1 H-NMR spectra for 10, 20, and 150 K. 
The intensities at different temperatures are not to scale; for 
details see text. 



more pronounced for the directions perpendicular to the 
Cu chain. The latter is in perfect agreement with the 
results obtained from the temperature dependence of the 
ESR linewidth and emphasizes the magnetic anisotropy 
in our system. 

The temperature dependence of the spin-lattice relax- 
ation rate, 1/Ti, is depicted in Fig. 12 for both probing 
nuclei and for the external magnetic field parallel to the a, 
b, and c axes. Here, we present our spin-lattice relaxation 
data in both the paramagnetic and the magnetically or- 
dered state below Tn. This way, a first microscopic proof 
for the 3D magnetically ordered state is given. Due to 
the very short spin-spin and also spin-lattice relaxation 
times of the order of 5 and 20 /us, respectively, in the 
temperature range close to Tn, our data are marked by 
large error bars in this temperature region and an accu- 
rate determination of Tn via 1/T. is difficult. From the 
overall temperature dependence of the spin-lattice relax- 
ation rate Tn was determined to be about (2.5±0.2) K at 
2T. 

While the longitudinal nuclear magnetization is well 
described by the standard expression of a nuclear spin / 



= 1/2 with a single T\ component for 207 Pb and for all 
field directions, for 'H an additional T\ component with 
a very short spin-lattice relaxation time of the order of ~ 
10 /xs needs to be taken into account for the whole tem- 
perature range for H || b. 55 Although it cannot be ruled 
out completely that the short Ti-component arises from 
either impurities in the sample or from a ^ background 
signal from e.g. teflon outside the coil, both its qual- 
itatively similar temperature dependence as well as its 
complete absence for the other two directions, strongly 
hint towards an intrinsic signal. However, the origin of 
this very fast- relaxing component for H \\ b is not clear 
up to now and needs further investigation. Thus, we 
solely concentrate on the single, longer T\ component in 
the following. 

Overall, as temperature is lowered, for both nuclei one 
finds a strong increase in l/Ti(T) with a sharp peak at 
Tm « 2.5(2) K, below which the spin-lattice relaxation 
rate decreases again as shown for all three principle di- 
rections (Fig. 12). This sharp peak is indicative of a 
transition into the 3D ordered state. The divergence of 
the relaxation rate l/T. due to the transition at ~ 2.5 K 
makes a further comparison to any low-dimensional (ID, 
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FIG. 12: (Color online) The spin-lattice relaxation rate, 1/Ti, 
of PbCuS04(OH)2 as function of temperature determined in 
an external field of 2 and 4T for 1 H and 207 Pb, respectively. 
While the 207 Pb-NMR signal was wiped out at low temper- 
atures due to very short spin-spin and spin-lattice relaxation 
rates, the 1 H-NMR signal could be observed in the whole tem- 
perature regime. Note that the average 1/T-value has been 
extracted for H || b due to overlapping 1 H-NMR lines. The 
lines are guides to the eyes. 
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2D) models in the paramagnetic regime difficult, as it 
masks the quasi- ID behavior of linarite at low tempera- 
tures even above Tn. However, it is crucial to emphasize 
that the 1 H-NMR investigations in the 3D ordered state 
have shown a sharp splitting into a multi-peak pattern, 
which is consistent with a non-trivial antiferromagnetic 
alignment of the Cu-spins below Tn. Correspondingly, 
the 1/Ti data below Tn, shown in Fig. 12 as closed cir- 
cles, reflect the average value of 1/Ti for all lines with 
the small distribution of values marked by error bars. 



VI. THEORETICAL ANALYSIS - ASSIGNMENT 
OF EXCHANGE INTEGRALS 




A. Saturation fields and Curie- Weiss temperature 

As a starting point, we assume at first that the inter- 
chain interaction is weak and that the inchain interaction 
(frustration) obeys the rigorous two-magnon bound state 
condition 14 a > a 3c = 0.3676776, i.e., it is within the 
quadrupolar region at the saturation field. Then, using 
the experimental value for the Curie- Weiss temperature 
CW « 0.5|Ji|(l — a) w 27 K, reported above, the es- 
timated ID saturation field of about H% ~ 5.5 T (JT C 3, 
see Fig. 3), and the ^-factor for H \\ a, i.e., g a = 2.34 
(see Sect. IV), we obtain for the ratio r of the former 
quantities for the a direction 



e, 



1-a 2 



9V 



tH.<, 



Aa 2 + 2a - 1 



3.122575. 



(6) 



Note, that the last equation is a rigorous relation valid 
in the commensurate 14,56 ID quadrupolar phase (QP) of 
the adopted spin model. Inverting Eq. (6), we obtain for 
the frustration ratio 



a(r) 



y/5 + 5/r + l/r 2 -! 
4+ 1/r 



0.36784 



(7) 



still in the quadrupolar phase, although very close to the 
border to the octupolar phase at a3 C . From CW we 
obtain J x = -85.4 K and J 2 = 31.4 K. 

In the opposite limit of strong enough antiferro- 
magnetic interchain coupling, where multipolar effects 
disappear, 39 the former can be described at least approx- 
imately by the well-known one-magnon theory in the case 
of separable in-chain and interchain contributions to the 
total 3D (2D) saturation field. 14 As a result one arrives 
at 



0.5 



2a- 1 + 0.125/a' 
0.5 



(8) 



a 



4r + 1 
0.3639, 



2r + 1 - 0.5p + v / (3-2p)r+(l-0.5p) 2 

(9) 



FIG. 13: (Color online) Maximum position T™ ax of the spin 
susceptibility x(T) measured in units of Ji as a function of 
1/a, i.e., viewing the J1-J2 model as an equivalent realization 
of two interpenetrating interacting AFM Heisenberg chains. 
In the inset the maximum position y™ ax of x(T) measured 
in units of |Ji| vs. the frustration parameter a according to 
DMRG calculations is shown (see text). The broken line de- 
scribes the asymptotic curve T™ ax = 0.641 J 2 = 0.641a|Ji|, 
well known from the Bethe-Ansatz-based solution for the un- 
frustrated AFM spin-1/2 Heisenberg chain. 57 



where /3 denotes the 2D interchain coupling in the basal 
plane measured in units of \Ji \ and p = JT 3D /i/ 1D — 1 rj 
0.3819 is given by the ratio of the 3D(2D) and the ID 
saturation fields, 7.6 T and 5.5 T, respectively, and the 
experimental value r = 3.122575 has been used. These 
numbers result in J\ = —86.7 K and J2 = 31.6 K. No- 
tice that the obtained value a is very close to the "two- 
magnon-value" estimated above and the value obtained 
from the analysis of the magnetization measurements (see 
Fig. 14). 



B. Spin susceptibility and magnetization curve 



Noteworthy, from the maximum position of the spin 
susceptibility at T™ ax = 4.9 ± 0.3 K in units of J 2 or 
|Ji| (Fig. 13) we can also estimate an a- value for our 
compound. Here, T™ ax is described with high precision 
by 



T~imax 
X 



(a) 



Jl\ 



E 

m— 1 



A,, 



1 

a-- 



with Ai = 
1.8487, A 5 



for 9/4>a>l/4, 



(10) 

-2.559, Ai = 



0.2778, A 2 = 1.7055, A 3 
= -0.6499, and A 6 = 0.0891. 
The expression given by Eq. (10) has been obtained 
from a fit of our TMRG-data for x{T) for strong and 
intermediate coupling. We note that Eq. (10) also de- 
scribes the experimental situation and paramctrizations 
suggested for Li2ZrCu04 with a = 0.3 34 as well as for 
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LiVCu0 4 with a = 0.75. 39 ' 58 In the opposite weak- 
coupling limit it is convenient to expand T™ ax around 
the limiting point of decoupled interpenetrating AFM 
Heisenberg chains in powers of 1/a = \ J\\/J 2 : 



X 



o) 



J-2 



= 0.641 



a™ ' 

rn—2 



for a > 1, 



(ii) 



-0.669, and D 5 = 



with D 2 = 0.0034, D 3 = 0.499, L> 4 
0.281. 

For linarite, Eq. (10) yields somewhat larger values for 
the same a, namely J\ = -107.5 K and J 2 — 39.5 K. We 
ascribe this difference to the presence of an antifcrro- 
magnctic interchain coupling ignored in the former esti- 
mation and an uncertainty in the determination of C w 
(Note that T™ ax is not affected by the interchain cou- 
pling within the random phase approximation adopted 
for its treatment at variance to the Curie- Weiss temper- 
ature Ocw) 

A completely different situation has been suggested 
by Y. Yasui et al. [22]. These authors arrived at a 
ymax ^ 3 g to 5 K somewhat lower than our values but 
at a negative C w ~ — 4 K, J\ = (— 13 ± 3) K, and 
an NNN exchange integral J 2 = (21 ± 5) K, yielding a 
frustration ratio a = 1.6 belonging to the weak-coupling 
region. Adopting their a-value and using their experi- 
mental T™ ax of about 5 K (see inset in Fig. 2 of Ref. 
22), one would arrive at J 2 « 8 K and at J\ « —5.2 K, 
which is inconsistent with (21±5) K and (—13 ± 3) K, 
respectively, claimed before. 

Next, we checked these preliminary values of the in- 
chain exchange integrals by additionally describing the 
low- and intermediate-field magnetization data at T = 
1.8 K for the a± direction (Fig. 14). The high-field re- 
gion has not been taken into account due to an addi- 
tional phase transition of yet unknown nature at low 
temperatures. The best fit at low fields is obtained for 
a ~ 0.365. Then, the corresponding Ji-value consistent 
with a ID saturation field of H s — 5.5 T becomes J\ = 
-89.5 K. These value are in accord with the values esti- 
mated above. Therefore they might be regarded as first 
realistic phcnomenological values despite the ignored but 
certainly present weak interchain and spin-anisotropy ef- 
fects, only briefly discussed below as an outlook for future 
studies. 

Finally, the analysis of the spin-susceptibility data 
nicely confirms the above results (Fig. 15). In the analy- 
sis the obtained experimental parameters for the (/-factor 
and the Curie- Weiss constant have been used. The best 
matching between experimental data and theoretical cal- 
culations results in magnetic exchange parameters J\ = 
-94 K and J x = -101.2 K (a = 0.36) for the a and b 
direction, respectively, again highlighting the magnetic 
anisotropy in our system. 

In summary, within the simplest isotropic spin chain 
model our values for J\ w —97 ± 10 K and J 2 w 36 ± 4 K 
exceed significantly those determined previously from fits 



of the T-dependent susceptibility in the high-T regime, 
where only Jj = -30 K (-13 K) and J 2 = 15 K (21 K), 
respectively, have been extracted. 22 ' 33 In addition, we 
also found a considerably smaller frustration parame- 
ter a « 0.37 instead of 0.5 and 1.62, respectively, 22 ' 33 
which puts linarite closer to the ID critical point at 
0.25 with consequences for a weaker critical antifcrromag- 
netic interchain coupling for ordered multipolar phases at 
T = 0. 18 The former might be masked by non- negligible 
impurity contributions in accord with a 3D analysis of 
the anisotropic susceptibility data (see below). 



C. Aspects of interchain coupling and symmetric 
exchange anisotropy for x(^) 

The interchain coupling (IC) has been taken into ac- 
count in the frame of the frequently used random phase 
approximation (RPA): 57 ' 59 



XlD 



l + fc X io(T) 



(12) 



where the one-dimensional susceptibility, Xm(T), has 
been calculated applying the TMRG-method, the inter- 
chain coupling k = (g/2) 2 J2Jic/\Ji\ yields a tempera- 
ture independent single parameter and the summation 
runs over all interchain couplings. The fits to the data 
shown in Fig. 15 result in slightly different IC param- 
eters k = —0.081 and -0.1 for the a and b direction, re- 
spectively, corresponding to a weak FM IC of a few K. 
In principle, a small underestimation of x(T) within its 
maximum region can be removed adopting also an impu- 
rity contribution described by a Curie or a Curie- Weiss 
law as frequently used in the literature for other quasi- 
1D compounds (see e.g. Refs. 57,60). Despite the uncon- 
trolled deviations introduced by the RPA, this way, the 
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FIG. 14: (Color online) Fit for the magnetization vs. external 
field of PbCuS0 4 (OH) 2 at T = 1.8K and for H || a±. 
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FIG. 15: (Color online) Analysis of the magnetic susceptibil- 
ity, x(T), within the isotropic J1-J2 model. The fitted solid 
lines result in slightly different IC parameters k = —0.081 to 
-0.1 (see text) corresponding to a weak FM IC of a few K. (a) 
H || a, (b) H || b axis. The g- factors are taken from Sect. IV. 



determination of the interchain coupling is not unique 
and a (sizable) impurity contribution might mask weak 
antifcrromagnetic interchain interactions. 

A weak AFM IC of the order of J lc w 4 K follows 
also from the positive difference between the estimated 
3D saturation fields and its ID counter part of « 6 T. 14 
Alternatively, the high-field magnetization should be af- 
fected significantly by changed exchange integrals in the 
new high- field phase observed above ~ 6 T. Inelastic neu- 
tron data such as in Ref. [61] might be useful to achieve 
a more precise assignment of the details such as the prin- 
ciple nature of the interchain coupling in our compound. 

Due to the increase of a possible, small impurity con- 
tributions with decreasing T in our compound, the max- 
imum positions of the impurity corrected susceptibilities 
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FIG. 16: (Color online) Magnetic susceptibility for a mag- 
netic field applied along the suggested easy-axis (b axis) fitted 
within the ID anisotropic J1-J2 model based on TMRG cal- 
culations supplemented by zero (upper panel) and finite anti- 
ferromagnetic (lower panel) interchain interactions treated in 
the RPA. The adopted easy-axis anisotropy is measured by 
the dimensionless parameter A > 1 (see text). 



are upshifted by a few tenths of a K. Naturally the ef- 
fect is largest for the 6-axis susceptibility (shift ~ 0.5 K). 
Nevertheless, a significant anisotropy remains also for 
these "corrected" susceptibilities which points together 
with the obtained largest exchange integral J\ for the b 
axis to a corresponding easy-axis assignment. To improve 
further the theoretical analysis, a detailed consideration 
of an anisotropic J\-3i model would be desirable which, 
however, is outside the scope of the present paper. For a 
first estimate see Fig. 16 and the short discussion given 
below. 
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Noteworthy, the quality of the fits can be improved 
without adopting a Curie- Weiss like impurity contri- 
bution, if instead a model with symmetric exchange 
anisotropy is applied. According to previous theoretical 
investigations of edge-shared cuprate chains with a fer- 
romagnetic nearest neighbor exchange coupling, an easy- 
axis exchange anisotropy is expected. 62,63 Therefore the 
first term in the isotropic spin-Hamiltonian should be re- 
placed by 

JiSiSi+i -> Ji (SfSf +1 + S*S* +1 ) + JiAiSfSf +1 , 

with Ai > 1 (see Fig. 16). Such an easy-axis exchange 
anisotropy affects significantly the low-temperature be- 
havior of x{T) and the saturation field in the low a- 
rcgion of interest. Following Rcf. 14, but ignoring weak 
anisotropy effects for the next-nearest neighbor exchange 
(i.e., setting A2 = 1), we rewrite the leading term in the 
expression for the ID saturation field valid in the two- 
magnon case at weak interchain coupling as 

S/i B # fl (Ai) = |Ji (2a - Ai + 0.5A 2 / [Ai + a]) , (13) 

and analogously in the one-magnon case: 

fl /i B #«(Ai) = I Ji| (2a — Ai + 0.125/a) . (14) 

The effect of the easy-axis anisotropy under considera- 
tion on the in-chain exchange coupling J\ is illustrated 
in Fig. 17 for the case we are interested in here, namely, 
when one extracts the J\ -value from a given experimental 
saturation field H s . 

In the presence of an antiferromagnetic interchain cou- 
pling the enhancement effect for the corresponding renor- 
malization of J\ is a bit less dramatic but nevertheless 
significant. Such a behavior is in accord with a similar 
effect found for the magnetic susceptibility data fitted 
by the isotropic and the aniosotropic model as shown in 
Figs. 15,16. Thus, within such a scenario the observed 
increase of | J\\ by more than 30 K becomes rather nat- 
ural. These results are in reasonable agreement with re- 
cent L(S)DA+[7 calculations, based on a refined crys- 
tallographic structure of PbCuS04(OH)2 including the 
corrected proton positions, 38 see below. Note that the 
range of validity of the expression given by Eq. (13) with 
respect to the formation of three-magnon bound states 
is under investigation at present. However, qualitatively 
the same behaviour is also expected for the octupolar 
three-magnon and other multipolar phases. A system- 
atic study of exchange anisotropy effects including also 
J2 as well as the interchain coupling is postponed for a 
future study. 



account the strong Coulomb repulsion, U, at the Cu site. 
Since we observed a sizable dependence of the resulting 
exchange parameters from the H position, 64 our calcula- 
tions are based on the recently refined crystal structure 
of Rcf. 38. For a screened Coulomb repulsion U = 7 eV, 
which is in the typical range of U values that were suc- 
cessfully applied to related Cu-0 systems, 36,61:65 and the 
usual value of the Hund's rule coupling J = 1.0 cV which 
enters the L(S)DA+{7 calculational scheme, we obtain a 
frustration ratio a ps 0.32. This value is in reasonable 
agreement with the value a ~ 0.36 derived from the ex- 
perimental data (compare Tab. 1), in view of possible 
renormalizations due to a non-negligible spin-lattice cou- 
pling in the non-adiabatic limit or intermediate case and 
of possible quantum-effect caused by the zero-point mo- 
tion of the light hydrogen nuclei ignored in all density 
functional approaches. 

The resulting frustration ratio a depends only mod- 
erately on U within the physically reasonable range be- 
tween 5 and 8 eV (Fig. 18). The respective calculated 
exchange integrals (U = 7 eV) are J\ = -133 K and J2 
= 42 K. These numbers are in good agreement with the 
isotropic parts obtained within the easy-axis model for 
the fit of the magnetic susceptibility with H \\ b reported 
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D. Density functional theory: L(S)DA + U 

To probe our new parameter set with respect to a mi- 
croscopic picture, we carried out DFT band structure cal- 
culations within the LSDA+f/ scheme which takes into 



FIG. 17: (Color online) Value of the isotropic part of the 
ferromagnetic nearest-neighbor exchange integral Ji in units 
of the corresponding value in the isotropic model vs. the 
dimensionless anisotropy factor Ai > 1 as derived from Eqs. 
(13) and (14). 
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above: — Ji = 138 K and J2 = 49. 7K. In other words, 
also the LSDA+J7 derived exchange integrals clearly sup- 
port a scenario with significantly larger J- values com- 
pared to those proposed in Refs. 22,33. Taking into ac- 
count the different interchain couplings Ji Cj o, Jic,i, and 
Jic,2 (see Fig. 1) as well as the respective number of 
neighbors we can estimate an effective interchain cou- 
pling Ji C ~ 7 K. This is also consistent with the analysis 
of x(T) within the easy-axis model which yields w 5 K. A 
detailed electronic-structure study of the title compound 
under consideration, in particular its structural stability 
and the sizable dependence of the coupling parameters 
from the details of the crystal structure, will be published 
elsewhere. 64 



VII. CONCLUSION 

In conclusion, we have performed a detailed combined 
experimental and theoretical study of PbCuS04(OH)2 in 
the paramagnetic regime in fields up to saturation. The 
saturation field, probed via magnetization studies at low 
temperatures, is anisotropic for the three principal crys- 
tallographic axes ranging between 7.6 and 10.5 T. ESR 
and NMR measurements further prove that the static 
susceptibility is dominated by the intrinsic spin suscepti- 
bility. The Knight shift as well as the broadening of the 
linewidth in the ESR and NMR data at elevated tem- 
peratures indicate a frustrated system with the onset of 
anisotropic magnetic correlations far above the magnetic 
ordering temperature, Tn = 2.75(5) K. Our experimen- 



axis 


a 


b 


9 (ESR) 


2.34 


2.1 


Ocw (K) 


27 


27 


H B (T) from 


5.5 




ID M(H,T= 1.8K) 






-Ji (K) (IHM, X (T)) 


94 


101 


J 2 (K), (IHM, X (T)) 


33.8 


36.4 


-Ji (K) (IHM, M(H)) 


89.5 




J 2 (K) (IHM, M(H)) 


32.7 




-Ji (K), (AHM) 




138 


a (IHM, M(H)) 


0.365 




a (AHM) 




0.36 


-Ji (K), (LSDA+U) 


133 


133 


J 2 (K), (LSDA+U) 


42 


42 



TABLE I: Calculated magnetic exchange interactions Ji and 
Ji along the chain using a ID approach based on the analy- 
sis of saturation fields at T = 0, low- and intermediate field 
magnetization data within the DMRG, and LSDA + U calcu- 
lations, together with the experimental values for the g-factor, 
H s , and Qcw of PbCuS04(OH)2; for details see text. The 
abbreviations IHM and AHM denote isotropic and anisotropic 
Heisenberg model, respectively. 



tal data are analyzed both in ID as well as quasi-lD 
approaches based on DMRG and hard-core boson cal- 
culations, yielding values for the exchange interactions 
J\ ~ —100 K and J2 ~ 36 K along the chain as well as a 
weak interchain coupling of a few K, leaving room for a 
quadrupolar phase (i.e., spin nematics) in experimentally 
accessible magnetic fields. 

In view of the small absolute value of the saturation 
field and the rich manifold of field-induced phases (partic- 
ularly along the chain direction), PbCuS04(OH)2 reveals 
itself as a promising material for additional investigations 
in the future. In particular, linarite possibly realizes a 
prototype for investigating the recently predicted spin 
multipolar order close to the saturation field, H s , of such 
frustrated "ferromagnetic" spin-chain compounds. 

From a theoretical point of view, such a systematic 
study of one compound on high-quality single crystals 
might be very helpful for a future unified theory of frus- 
trated edge-shared cuprates in order to understand what 
is really puzzling in each compound of that increasing 
family and what some of these compounds have in com- 
mon. In particular, the question which of their proper- 
ties can be used best for the determination of a specific 
model parameter from experimental data seems to be im- 
portant. 
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FIG. 18: Frustration parameter vs. U from L(S)DA+(7 for 
PbCuS04(OH)2. The Hund's rule exchange at Cu sites has 
been fixed at J = 1.0 eV for all calculations. 
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